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Elliptic-Cylindrical Wavelets: The Mathieu Wavelets 
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Abstract 

This note introduces a new family of wavelets and a multiresolution analysis, which exploits 
the relationship between analysing filters and Floquet’s solution of Mathieu differential equa¬ 
tions. The transfer function of both the detail and the smoothing filter is related to the solution 
of a Mathieu equation of odd characteristic exponent. The number of notches of these filters 
can be easily designed. Wavelets derived by this method have potential application in the fields 
of Optics and Electromagnetism. 
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1 Introduction 

In 1868, the French mathematician E. Leonard Mathieu introduced a family of differential equations 
nowadays termed Mathieu equations in his “memoir on vibrations of an elliptic membrane^' [T]. 
Mathieu’s equation is related to the wave equation for the elliptic cylinder. Mathieu is notably 
remembered for his discovery of sporadic simple groups [2]. This paper is particularly concerned 
with the canonical form of the Mathieu Equation. Eor o € M, G C, the Mathieu Equation is given 
by ^ 

+ {a -2q cos{2uj))y = 0. (1) 

The Mathieu equation is a linear second-order differential equation with periodic coefficients. This 
equation was shown later to be also related to quantum mechanicals; the parameters a and q denote 
the energy level and an intensity, respectively. Eor g = 0 it reduces to the well-known harmonic 
oscillator, a being the square of the frequency [3]. The solution of ([T]) is the elliptic-cylindrical 
harmonic, known as Mathieu functions. In addition to being theoretically fascinating, Mathieu 
functions are applicable to a wide variety of physical phenomena, e.g., diffraction, amplitude dis¬ 
tortion, inverted pendulum, stability of a floating body, radio frequency quadrupole, and vibration 
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in a medium with modulated density [3]. They have also long been applied on a broad scope of 
waveguide problems involving elliptical geometry, including: (i) analysis for weak guiding for step 
index elliptical core optical fibres [5], (ii) power transport of elliptical waveguides PE!, (iii) evalu¬ 
ating radiated waves of elliptical horn antennas [S] , (iv) elliptical annular microstrip antennas with 
arbitrary eccentricity [9|, and (v) scattering by a coated strip m- 

The aim of this paper is to propose a new family of wavelets based on Mathieu differential 
equations. Wavelets are a well-known tool for differential equation solving [nHE]. However, in 
this work, we show another connection between wavelets and differential equations: the design of 
new wavelets from the solution of a differential equation. 

2 Mathieu Equations 

In general, the solutions of ([U are not periodic. However, for a given q, periodic solutions exist for 
infinitely many special values (eigenvalues) of a. For several physically relevant solutions y must 
be periodic of period vr or 27 r. It is also convenient to distinguish even and odd periodic solutions, 
which are termed Mathieu functions of first kind. One of four simpler types can be considered: 
Periodic solution (vr or 27 r) symmetry (even or odd). For q ^ 0, the only periodic solution y 
corresponding to any characteristic value a = ar{q) or a = br{q) has the following notation: 

Even periodic solution 


cer{co,q) = Ar^m cos mcj for a = ar{q), ( 2 a) 

m 

Odd periodic solution 

ser{io,q) = sin mo; for a = 6 ^( 9 ), ( 2 b) 

m 

where the sums are taken over even (respectively odd) values of m if the period of y is vr (respectively 
2vr). Given r, we denote henceforth by Am, for short. Elliptic cosine and elliptic sine functions 
are represented by ce and se, respectively. Interesting relationships are found when g —)• 0, r / 0 [T3] : 

limcer(a;,g) = cos(ra;), limser(a;,g) = sin(ra;). (3) 

One of the most powerful results of Mathieu’s functions is the Floquet’s Theorem m- It states 
that periodic solutions of ([T]) for any pair (a, q) can be expressed in the form 

y{uj) = Fy{uj) = e^^'^P{u) or 
y{u) = F,{-u) = e-^"‘^P(-w), 

where w is a constant depending on a and q and P{-) is vr-periodic in a;. The constant u is called the 
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characteristic exponent. If v is an integer, then F^{uj) and Fy{—uj) are linear dependent solutions. 
Furthermore, y{u} + kn) = or y{oj + fevr) = e~^'^^'^y{oj), for the solution Fy{Ljj) or Fy{—uj), 

respectively. We assume that the pair (a, g) is such that | cosh(jV7r)| < 1 so that the solution y{oj) 
is bounded on the real axis M- The general solution of Mathieu’s equation (g G M, non-integer) 
has the form 

y(a;) = P{uj) + C2e-^"‘"P(-u;), (5) 

where ci and C 2 are arbitrary constants. 

All bounded solutions —those of fractional as well as integral order— are described by an 
infinite series of harmonic oscillations whose amplitudes decrease with increasing frequency. In the 
wavelet framework we are basically concerned with even solutions of period 2tt. In such cases there 
exist recurrence relations among the coefficients |14j : 

(a - I - g)Ai - gAs = 0, 

( 6 ) 

(a - rr?)Am - q{Am -2 + ^m+ 2 ) = 0, m > 3, m odd. 

In the sequel, wavelets are denoted by and scaling functions by </>(t), with corresponding 
spectra 'I'(a;) and <I’(a;), respectively. 

3 Mathieu Wavelets 

Wavelet analysis has matured rapidly over the past years and has been proved to be invaluable 
for scientists and engineers m- Wavelet transforms have lately gained extensive applications in 
an amazing number of areas [TS]. The equation (/i(t) = ~ ^)) which is known as 

the dilation or refinement equation, is the chief relation determining a Multiresolution Analysis 
(MRA) dS]. 

3.1 Two Scale Relation of Scaling Function and Wavelet 

Defining the spectrum of the smoothing filter {hk\ by H{u) = central 

equations (in the frequency domain) of a Multiresolution analysis are |20j : 

^>( 0 ;) = R (I) $ (I) and T(a;) = G (^1) (^1) , (7) 

where G{aj) = , is the transfer function of the detail filter. 
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The orthogonality condition corresponds to |20j : 


F1(0) = 1 and H{tt) = 0, 

(8a) 

|R(a;)|2 + |R(u; + 7r)|2 = l, 

(8b) 

H{u) = -e-^^G*{uj + TT). 

(8c) 


3.2 Filters of a Mathieu MRA 

The subtle liaison between Mathieu’s theory and wavelets was found by observing that the classical 
relationship 

T(u;) = e-^W2^*(|-7r)c^(|) (9) 

presents a remarkable similarity to a Floquet’s solution of a Mathieu’s equation, since H{u}) is a 
periodic function. 

As a first attempt, the relationship between the wavelet spectrum and the scaling function was 
put in the form: 

i- - it) ■ (10) 

Here, on the second member, neither v is an integer nor H{-) has a period tt. By an appropriate 
scaling of this equation, we can rewrite it as 

- tt). (11) 

Defining a new function = T(4a;)/<h(2a;), we recognise that it has a nice interpretation in 

the wavelet framework. First, we recall that 41 ( 0 ;) = G (f)‘h (^) so that ^{2uj) = 

Therefore the function related to Mathieu’s equation is exactly Y{ui) = G{2uj). Introducing a new 
variable z, which is defined according to 2z = 2u — vr, it follows that —Y (^ + f) = e~^‘^^H*{2z). 
The characteristic exponent can be adjusted to a particular value u, 

_ ^-j(u-2)zY = e-^''^H*{2z). (12) 

Defining now P{-z) = H*{2z) = Zlfcez where C 2 k = we figure out that the right- 

side of the above equation represents a Floquet’s solution of some differential Mathieu equation. 
The function P{-) is vr-periodic verifying the initial condition F’(O) = = 1, as expected. 

The filter coefficients are all assumed to be real. Therefore, there exist a set of parameters (og, qg) 
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such that the auxiliary function 


is a solution of the following Mathieu equation: 


(13) 


dz^ 


+ (OG - 2qG cos(22;))yi, = 0, 


(14) 


subject to yjy(O) = —y(7r/2) = —G{tt) = —1 and cos(7rz^) — y^('/r) = 0, that is, yi/(vr) = (—1)^. 

In order to investigate a suitable solution of (jl4p . boundary conditions are established for pre¬ 
determined a, q. It turns out that when v is zero or an integer, a belongs to the set of characteristic 
values ar{q). Furthermore, v = r \s associated with ar{q). The even (27r-periodic) solution of such 
an equation is given by: 

ceu{z,q) 


yu{z) = -- 


(15) 


cei.(0,g) 

The Yjy[uj) function associated to yv{z) and related to the detail filter of a “Mathieu MRA” is thus: 


Y,{u) = G,(2w) = 

ceu{d,q) 


(16) 


Finally, the transfer function of the detail filter of a Mathieu wavelet is 

— oiT-2 


CQv{d,q) 


(17) 


The characteristic exponent iz should be chosen so as to guarantee suitable initial conditions, i.e., 
Gi/(0) = 0 and Gu^tt) = 1, which are compatible with wavelet filter requirements. Therefore, iz must 
be odd. It is interesting to remark that the magnitude of the above transfer function corresponds 
exactly to the modulus of a elliptic sine m- 


\G^{uj)\= se^ (^^,-q'^ /ceu{0,q) 


(18) 


The solution for the smoothing filter H{-) can be found out via QMF conditions [19], yielding: 


U C6;< 

ce,,(0,g) 

In this case, we find Hy{'K) = 0 and 

\Hy{oj)\ = / ceu{d,q) 


(19) 


( 20 ) 


Given q, the even first-kind Mathieu function with characteristic exponent v is given by cey{uj, q) = 
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Figure 1: Magnitude of the transfer function for Mathieu multiresolution analysis filters: smoothing 
filter \Hy{uj)\ (solid line) and detail filter \Gu{uj)\ (bold line) for a few Mathieu parameters, (a) 
y = 2,, q = 2>, a = 9.915506290452134; (b) = 5, g = 15, a = 31.957821252172874. 


^ 2 Z+i cos(2/ + l)w, in which cei/(0,g) = X]/^o^ 2 Z+i- The G and H filter coefficients of a 
Mathieu MRA can be expressed in terms of the values {A 2 i+i}i£z of the Mathieu function as: 


hf _ T|a-^|/2 fff y 7l|2;+^_2|/2 

V2 ce,(0,g) ^2 ^ ^ ce,(0,g) 


( 21 ) 


It is straightforward to show that VZ > 0. The normalising conditions are 

71 Er=-oc K = -1 and EZ-ooi-^)’^K = 0. 


4 Examples 

Illustrative examples of filter transfer functions for a Mathieu MRA are shown in Fig. [U for v = 3 
and 5, and a particular value of q (numerical solution obtained by 5-order Runge-Kutta method). 
The value of a is adjusted to an eigenvalue in each case, leading to a periodic solution. Such 
solutions present a number of v zeroes in the interval |w| < vr. We observe lowpass behaviour 
(for the filter H) and highpass behaviour (for the filter G), as expected. Mathieu wavelets can be 
derived from the lowpass reconstruction filter by the cascade algorithm. Infinite Impulse Response 
filters (HR) should be applied since Mathieu wavelet has no compact support. However a Finite 
Impulse Response (FIR) approximation can be generated by discarding negligible filter coefficients, 
say less than 10“^*^. In Fig. [SJ emerging pattern that progressively looks like the wavelet shape is 
shown for some couple of parameter a and q. Waveforms were derived using the Matlab wavelet 
toolbox. As with many wavelets there is no nice analytical formula for describing Mathieu wavelets. 


5 Conclusions 

A new and wide family of elliptic-cylindrical wavelets was introduced. It was shown that the transfer 
functions of the corresponding multiresolution filters are related to Mathieu equation solutions. The 
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(a) 







(b) 

Figure 2; FIR-Based Approximation of Mathieu Wavelets as the number of iteration increases (2, 
4, and 6 iterations, respectively). Filter coefficients holding \h\ < 10“^® were thrown away (19 
retained coefficients per filter in both cases), (a) Mathieu Wavelet with u = 3 and q = 3 and (b) 
Mathieu Wavelet with u = 5 and q = 15. 
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magnitude of the detail and smoothing filters corresponds to first-kind Mathieu functions with odd 
characteristic exponent. The number of zeroes of the highpass |G(w)| and lowpass \H{u!)\ filters 
within the interval |a;| < vr can be appropriately designed by choosing the characteristic exponent. 
This seems to be the first connection found between Mathieu equations and wavelet theory. It opens 
new perspectives on linking wavelets and solutions of other differential equations (e.g. Associated 
Legendre functions). 

Although there exist plenty of potential applications for Mathieu Wavelets, none are presented: 
we just disseminate the major ideas, letting further research to be investigated. For instance, this 
new family of wavelets could be an interesting tool for analysing optical fibres due to its “elliptical” 
symmetry. They could as well be beneficial when examining molecular dynamics of charged particles 
in electromagnetic traps such as Paul trap or the mirror trap for neutral particles |21[|22] . 
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